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Abstract 

A two-components isospin lattice gas model has been employed to study 
the liquid-gas phase transition for asymmetric nuclear matter. An additional 
degree of freedom, namely, the asymmetry parameter a has been considered 
carefully for studying the phase transition. We have shown that under the 
mean field approximation, the liquid-gas phase transition given by this model 
is of second order. The entropy continues at the phase transition point. The 
binodal surface is addressed. 
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An one-component lattice gas model with one type of atom has been found to be success- 
ful for describing the liquid-gas phase transition (LGPT) [0. An isospin lattice gas model 
which assumes a lattice with each site either being vacant or occupied by proton and neu- 
tron has been employed by many authors to study the LGPT of symmetric nuclear matter 
recently [2-6]. This model is mapped into a spin-1 Ising model. It has been proven that the 
LGPT of this model is of first order under the Bragg-Williams mean field approximation. 
The chemical potential continues at the phase transition point but its first order deriva- 
tives, namely, entropy and volume, are discontinuous. This result is reasonable if we notice 
that the number of protons equals to the number of neutrons for symmetric nuclear 
matter and the symmetric isospin lattice gas model is, in fact, a system with one-component. 

However, as was pointed out by ref.[7-12], the LGPT for a multi-components system 
is of second order, i.e. the entropy continues and the second order derivatives of chemical 
potential are discontinuous. This is because the degrees of freedom increase according to the 
Gibbs phase rule, and then the binodal surface has greater dimensionality. Since the space 
occupied by nuclear matter addressed by ref.[7-12] is continuous, the above conclusion is of 
course limited to a continuous manifold of space structure. It is of interest to extend this 
study to a two-components system with lattice space structure. 

The isospin lattice gas model for asymmetric nuclear matter suggested by ref.[13] is a 
typical two-components model because 7^ A^_, the number of protons and neutrons are 
independent. We will employ this model to reconsider this problem. In conflict with ref. [13], 
we will show the LGPT is of second order for two-components isospin lattice gas model and 
the binodal surface will end at two points with a = and a = 1 where a is the asymmetric 
parameter 



At these two end-points the LGPT becomes first order because the nuclear matter reduces 
to symmetric and neutron matter, and the lattice gas model reduces to one-component. 
The Hamiltonian of lattice gas model is[13] 



where < ij > denotes a nearest neighbor pair, h is a constant and the interaction strength 
parameter Jij satisfy 
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The isospin is 




— 1 for neutron 
for vacancy 
+1 for proton 
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The Hamiltonian can be written as 
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where A^++, and A^_| represent the nearest neighbor pairs of proton-proton, neutron- 
neutron and proton-neutron respectively, and the interaction strength parameter Jg for pp 
and nn pairs, and Jd for pn pairs. Introducing 
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where r denotes the relative emptiness and s plays the similar role as that of asymmetric 
parameter a. s is not a new parameter. We introduce it for comparing with ref. [13] only. 
By using the Bragg- Williams mean field approximation 
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where 7 denotes the number of nearest neighbors of any given site, and ^"yN is the total 
number of pairs. We can rewrite Eq.(|^) as 



where 
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The grand partition function of our system is 

r,s r,s 

where /3 = the factor 
g{r,s,N) 
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the fugacity z+ = e^^+ and Z- = e^^~ for proton and neutron respectively. For two- 
components system, in general, the chemical potential of proton /i+ does not equal to the 
chemical potential of neutron /i_ because 7^ A^_. According to the Gibbs phase rule, the 
degrees of freedom of a system is 



/ = k-ip + 2 



(12) 



where k and cp are the number of components and phases respectively. A two-components 
system has one more degree of freedom than that of one-component. We choose the asym- 
metric parameter a (or s) as the additional degree of freedom for two-components lattice 
gas model as that of nuclear matter models with continuous space structure [7-12]. There- 
fore, we have two independent parameters, namely r (corresponding to the density) and s 
(corresponding to the asymmetry of protons and neutrons) in a two- components lattice gas 



3 



model for asymmetric nuclear matter. This is the basical difference between our calculation 
and that of ref. [13]. 

In the thermodynamic limit, the sum in Eq.(|l0|) can be replaced by its most dominant 
term [3,13,14]. Since we have two independent parameters r and s, the extreme values (f, 
s) satisfy 

d In Qg 

\r=f,s=s = (13) 

and 



ainQc, 
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(14) 



We must use two equations to determin the extreme values of r and s. Using Stirling's 
formula and Eqs.(|T3|) (|1^, we obtian 

(3^+ = log z+ = - log f + log ^ - 2PC2 (1 - r) - 2pCis - I3h (15) 



= log 2_ = - log f + log - — ^ — - - 2(3C2 (1 - f) + 2/5Cis + (3h (16) 

Only the maximum of the extreme values gives us the most dominant contribution of Eq.(|lO|) 
and corresponds to stable state. Taking /i = to neglect the effect of the "external field", 
we find Eq.(|T5|) reduces to Eq. (p!6[) and /U+ = /i_ when s = 0, the asymmetric nuclear mat- 
ter becomes symmetric one and two-components becomes one- component. The symmetric 
nuclear matter is the specific case of our theory. 

By means of the grand partition function, we can obtain other thermodynamical quan- 
tities such as pressure, baryon density and entropy per baryon easily. They are 

7-, /- - o\ -2 , M -\2 {l-f + s) (1 - f + s) (1 - f - s) (1 - f - s) 
Pgas (r, S, (3) = Cis" + C2(l - vf - ^ — log ^ ^ - ^ — ^ log ^ ^ 

^1 7 (1 — f-hs), (1— f— s), 

--\ogf + hs + ^ — '-\ogz+ + ^ — ^log2_ (17) 

p=l-f (18) 
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^ ^log^ - + - ^log^ ^ + flogf (19) 



Now we are in the position to discuss the LGPT of our model. The two-phase coexistence 
equations read 

4 = ^+ (20) 
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= (21) 



= (22) 

where subsripts of one phase L and V stand for hquid and gas respectively. The parameters 
Ci, C2 are determined as follows: to have an attractive nearest neighbor interaction, we 
have Js > and Jd < 0. In this case C2 is positive. We choose the values of Ci and C2 to 
reproduce the binding energy {—16MeV) for symmetric nuclear matter at saturation point 
when T = OK, and find 

Ci = 3.2MeV, C2 = IQ.OMeV 

Our numerical results are shown in figures(l-4). In Fig.l, we show the isotherms for fixed 
a = 0.50 but different temperature T = 6.00, 6.90, 7.95, 8.58 and 8.70MeV respectively. 
The stable states of isotherms corresponding to the most dominant term in the summation of 



Eq.(|TO]) is drawn by the solid curves. The dotted curve which encloses the transition region 
represents the slice of binodal surface with constant asymmetric parameter. The critical 
temperature for a = 0.50 asymmetric nuclear matter Tc = 8.58Mey. 

The chemical isobar of proton and neutron vs. a curves at fixed temperature T = 
7.50MeV and p = 1.12MeV fm~^ are shown in Fig.2. In this figure, all metastable and 
unstable states are neglected. The Gibbs conditions Eqs. (20-22) for phase transition demand 
equal pressures and equal chemical potentials for two phase with different a. The collection 
of all such «! {T,p) and a2 {T,p) form the binodal surface. The two desired solutions form 
the edges of a rectangle and can be found by using the geometrical construction shown in 
Fig.2. We find the end points of the phase boundary are just mapped at the edges of the 
rectangle. This situation is the same as that of the models with continuous space manifold 
[7-12]. 

The section of binodal surface at T = 7.50MeV are shown in Fig.3. We find that the 
end-points of binodal surface are fixed at a = and a = 1 respectively. The pressure 
and chemical potential unchange during phase transition. This is just the behavior of one- 
component system. For cases < a < 1, the physical behavior of Fig.3 can be explained 
as follows: assume that the system is initially prepared with a = 0.5 (gas phase), during 
isothermal compression, the two phase coexistence region is encountered at point A, and the 
liquid phase emerges at point B. The isothermal compression process continues to the point 
C of liquid-gas coexistence region, the density p+ and /9_ for protons and neutrons satisfy 

p_ = Xpl + (1 _ A) (0 < A < 1) (23) 

where A refers to the volume fraction of gas phase, when the process evolves to the point F, 
the gas phase disppears and the system becoms liquid. 

Finally, we investigate the entropy and the order of LGPT. The entropy per baryon as 
a function of pressure for different asymmetric parameter a are shown in Fig. 4. The states 
denoted by dotted curves in the intermedial states during the phase transition is evaluated 
according to 
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(24) 



where p = + p-, p"^ = p^ + p^, p^ = p^ + p^ and the conservation law Eqs. (|23|) . The 
curves in Fig. 4 show the order of phase transition transparently. For symmetric nuclear 
matter {a = 0) and for neutron matter {a = 1), the entropy discontinues at the phase- 
transitional point. The LGPT is of first order for one-component system. For asymmetric 
nuclear matter, < a < 1, the entropy continues at the phase-transitional point and the 
LGPT is of second order. 

In summary, contrary to ref.[13], an additional degree of freedom a which makes the 
chemical potential of proton does not equal to that of neutron has been considered for 
asymmetric isospin lattice gas model. We have shown that under Bragg- Williams mean 
field approximation, the LGPT in a two- components isospin lattice gas model for asymmetric 
nuclear matter is of second order. At the phase transition point, the entropy continues. The 
binodal surface. The binodal surface for fixed temperature T = 7.50MeV has been given. 
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FIGURE CAPTIONS 

Fig.l The isotherms for fixed a — 0.50 and different temperatures T — 6.00, 6.90, 7.95, 8.58 
and 8.70MeV respectively. 

Fig.2 Chemical potentials for proton and neutron as a function oi a a,t p — 1.12MeV and 
T = 7.50MeV. 

Fig.3 The section of binodal surface for T = 7.50MeV. 

Fig.4 Entropy per baryon vs. pressure for asymmetric parameter a = 0.0, a = 0.6, a = 0.8 
and a = 1.0 at T = 7.5MeV. 
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